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ABSTRACT 

We show that, for a low-mass planet that orbits its host star within a few tenths of an AU (like 
the majority of the Kepler planets), the atmosphere it was able to accumulate while embedded in the 
proto-planetary disk may not survive unscathed after the disk disperses. This gas envelope, if more 
massive than a few percent of the core (with a mass below 10M®), has a cooling time that is much 
longer than the time-scale on which the planet exits the disk. As such, it could not have contracted 
significantly from its original size, of order the Bondi radius. So a newly exposed proto-planet would be 
losing mass via a Parker wind that is catalyzed by the stellar continuum radiation. This represents an 
intermediate stage of mass-loss, occurring soon after the disc has dispersed, but before the EUV/X-ray 
driven photoevaporation becomes relevant. The surface mass-loss induces a mass movement within 
the envelope that advects internal heat outward. As a result, the planet atmosphere rapidly cools 
down and contracts, until it has reached a radius of order 0.1 Bondi radius, at which time the mass- 
loss effectively shuts down. Within a million years after the disk disperses, we find a planet that has 
only about ten percent of its original envelope, and a Kelvin-Helmholtz time that is much longer 
than its actual age. We suggest that this “boil-off” process may be partially responsible for the lack 
of planets above a radius of 2.51?® in the Kepler data, provided planet formation results in initial 
envelope masses of tens of percent. 

Subject headings: planets and satellites: composition, planets and satellites: formation, protoplanetary 
disks, planet - disk interactions 


1. INTRODUCTION 

The Kepler observatory has revealed the presence of 


iBorucki et al.112011 

Batalha et al.||2013 

Mullally et al. 

GO 

i-H 

o 

CM 

and low-mass 

5 20-Mfn) exoplanet 

s (IVlarcv et al. 

2014 

Dumusque et al. 2014 

Dressing et al. 2015). in tact 


this type of planet is so frequent mo st stars are thought 


to contain one (e.g. 

Fressin et al. 

2013 Petigura et al. 

2013 Silburt et al. 

20150. We nick-name these plan- 


mode of planet formation, allowing us to calibrate our 
planet formation models and finally understand this il¬ 
lusive problem. 

Unlike the terrestrial planets in our own solar system, 
the inferred masses and radii of the Kepler planets indi- 


(IWu & Lithwickl 

2013; Hadden & Lithwick 

2014 Weiss 

& Marcy' 2014 

tiogers| 2U14p. Even among ones that 


dynamical escape ma y have occurred and removed their 
primordia l envelopes (Owen fe Wu||2013 Lopez & Fort¬ 


ney 


2013). This indicates that, in most cases, planet 
assembly must have finished by the time gas-rich proto¬ 
planetary disks disperse after a few Myr. 

The size distribution of Kepler planets exhibits a puz¬ 
zling feature: there appears to be a concentration of plan¬ 
ets with sizes ~ 2.5i?®, with a steep fall-off in number at 
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smaller sizes (Petigura et al. 

2013 

Foreman-Mackey et 

al. 

2014 

Silburt et al. 

2U15j), although the quantitative 


a size of 2.5J?m corresponds to an envelope mass fraction 
of ~ 1% ( Wu fe Lithwick|2013| Wolfgang &; Lopez|2015 ). 


Why does nature favour such an envelope mass, when the 
environments (e.g., gas density, temperature) for planet 
formation may be highly diverse and time varying? 

There are currently two schools of thought for where 
these planets formed: either they are formed “in-situ”, 
close to their current small or bital separations (e.g. 


Hansen & Murray 2012 Chiang & Laughlin 2013), or 


they formed at larger sep arations and are then migrated 


to current positions (e.g. Ida & Lin 

2005 

2010 Boden- 

heimer & Lissauer 2014 

Raymond & Cossou 2014j). Even 


gas accretion may have occurred near the host stars. Ke¬ 
pler planets that are in high multiple systems appear to 
be so closely packed (e.g., the Kep ler-11 system, Lissauer| 


et al.||20fT Mahajan fc Wu|20l4 1, they would have been 


skirting dynamical instabilities in the final assembly had 
they no t been subject to substantial eccentricity dissi¬ 
pation (Pu &_Wu 2015 1. Their current low values of 
eccen tricities (|Wu fe Lithwick 2013| |Hadden & Lithwick 
20141 also attest to this dissipation! Gas is the most 
obvious source of this damping. 

Could the gas accretion process, at a distance of ~ 0.1 
AU from the host star, naturally give rise to the above- 
mentioned 1% envelope mass? The answer is no, at least 
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not according to cu r rent theories. In the theory of core¬ 
accretion (e.g. Pollack et al. 1996|), a planetary embryo 
embedded in a gas disk will quickly acquire a hydrostatic 
envelope extending out to its Bondi splierfQ with a ra- 

diulO 

(1) 

where M p is the proto-planet’s mass and c s is the isother¬ 
mal sound speed of the surrounding gas. The atmo¬ 
sphere, confined by an external disk pressure, has a mass 
that is roughly the background gas d ensit y times the vol¬ 
ume of the Bondi sphere (Rafikov 2006). If this is the 
final atmosphere that we observe, its mass should de¬ 
pend on gas density, temperature and embryo mass. The 
observed clustering around 1% mass fraction therefore 
comes as a surprise. 

Multiple works have suggested that the abov e mass es¬ 
timates ma y not be the final envelope mass. Ikoma & 
Hori (2012) argued that, as the gas disk gradually dissi- 


pates away, lifting of the pressure confi nement will slowly 
erode away the envelope. In contrast, Lee et al. (20141 
asserted that due to a lack of planetesimal bombardment 
and associated heating for embryos in the inner disk, the 
above mentioned hydrostatic envelope will cool and con¬ 
tract, allowing the embryo to accrete ever more gas. In 
fact, they found that even a 5planet may accrete 
so m uch gas as to undergo unstable run-away growth 
(Mizuno et al. 19781. We return to comment on these 
works later in the discussion. Instead, in this work, we fo¬ 
cus on a different aspect of the problem, namely, the case 
where the nascent disk disappears qui ckly (as indicated 


observationally, e.g. Koepferl et al. 2013), leaving the 


newly formed planet exposed to stellar irradiation and a 
vacuum boundary condition. We argue that this brings 
about vigorous mass-loss and rapid cooling, strongly im¬ 
pacting the final envelope mass. This process, happening 
in the final stage of low-mass planet formation, may po¬ 
tentially help explain the 1% convergence one observes 
in Kepler planets. 

This new stage of mass-lo ss is distinctive ly different 


from the process discribed in Ikoma & Hori (2012) and 


ration ( 

Lopez & Fortney 

2013 

Owen & Wu 

2013 

). As 

our calculations show (§3), it occurs because t 

he planet is 


highly inflated and because it is strongly illuminated by 
stellar photons. Both these conditions may be satisfied 
for low-mass planets that newly emerge out of the proto¬ 
planetary disk, due to their long thermal timescales (§2). 


T he m ass-loss takes the form of a Parker wind (Parker 
1958) that is heated by continuum stellar radiation (as 


opposed to only high energy photons). But the ultimate 
energy source is the internal heat of the planet. As a re¬ 
sult of envelope removal, the planet cools and contracts 
quickly (§4). We consider the implications of this work 
in §4 and conclude in §5. 

2 In the inner disk, the Bondi radius is smaller than the Hill 
Spher e for the core masses we are concerned about (e.g. Rafikov] 

|20M|. 1 - 1 

It is yet unclear exactly where the proto-planet’s envelope 
ends and where the disk be gins. Back ground shear in the disk can 
affect the transition rad ius (Lissaucr et al. 2009]l, as can 3-D flow 

02015, Fung et af. 20151. 


dynamics (Ormel et al. J201 
complications nere. 


Fung et al.||2uii3|fT We ignore these 


In this work we concentrate purely on the case of a 
rocky core surrounded by a H/He envelope, i.e. the struc¬ 
ture most likely to arise from planets forming close to 
their star in a gaseous disc, although we note many other 
structures are s till consistent with the observed planet 
population (e.g. Rogers & Seagcr 2010). 

2. INITIAL CONDITIONS AND PRESENCE OF MASS-LOSS 

Planets can lose their atmospheres over time. This 
process is more extreme when the planets are strongly 
irradiated and highly inflated. We argue here that this 
may be just the thermal state the Kepler planets found 
themselves in, when their parent disks dissipated. The 
process we are interested in occurs within 10 4 — 10 5 yrs 
after the disk dispersal, before the EUV/X-ray photoe¬ 
vaporation has had much of an eff ect. T he latter occurs 
on a time-scale of 10 8 yrs ( Owen &; Wu||2013 ). 

2.1. A wind is driven 

The black-body temperature at a distance a from a 
star with surface temperature T* and radius 1?* is 


T eq = 886 K 


T* 


5800 K 


R * 


1 R f 


J vo.i 


AU 


- 1/2 


( 2 ) 

This is also the photospheric temperature of a planet 
when its internal luminosity is subordinate to stellar irra¬ 
diation and when it is exposed to direct sun-light. In this 
paper, we assume that the gas above the photosphere, in 
particular, the extended atmosphere that is the outflow 
we are investigating, will also be heated to roughly this 
temperature. This could happen due to the absorption 
of continuum photons by the dust grains embedded in 
the outflow, when they are present, or by absorption of 
ionizing photons by the gas molecules when dust is ab¬ 
sent. The planet distances of concern lie beyond the 
dust sublimation radius. So we expect dust formation to 
be efficient in the outflow. In this work, we do not at¬ 
tempt to model the temperature profile in the extended 
atmosphere, a short-coming we hope to remedy in future 
works. We do remark that the condition of isothermal 
wind sets an upper bound to the mass-loss rate: the en¬ 
ergy advected by the wind should not exceed stellar en¬ 
ergy deposit in the wind, or else the isothermal condition 
will break. 

An isothermal atmosphere has an outwardly increas¬ 
ing scale height and can only be in hydrostatic equi¬ 
librium if the ambient pressure is sufficiently strong to 
confine it. This confining pressure was provided by the 
proto-planetary disk while it was around, but once it 
dissipated, the only candidate is the stellar wind, which, 


(e-g- 

Wood et al. 

2002 

), falls far too short to be of use. 

So the isotherma 

atmosphere is not hydrostatic but in- 

stead harbours a so-called “Parker wind” ( 

Parker] 1958) 


with a sonic radius at: 


n _ gm p 

R s 2c? Rb ' 


(3) 


Note the fact that the sonic radius and the Bondi ra¬ 
dius are identical is not accidental, but represents the 
symmetry between outflow and accretion. 

Energetically speaking, the stellar continuum radia¬ 
tion, which maintains the isothermality of the wind, is 
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Fig. 1.— The internal luminosity of a planet with a 5M® rock 
core and a surface temperature of 900 K (roughly the blackbody 
temperature at 0.1 AU) as a function of the envelope mass fraction. 
The colour of the lines indicates the radius of the planet in units 
of the Bondi-radius . The three horizontal lines lab el the expected 
accretion luminosities from solids in a MMSN disk (|Rafikov|2006|). 
At these values, solid accretion easily overwhelm the cooling con- 
traction of the proto-planets. 


the reason for the mass-loss. However, as we will dis¬ 
cover in later parts, the mass-loss brings about a rapid 
gravitational contraction of the planet. The binding en¬ 
ergy released during this contraction fuels the mass-loss. 
There are thus two energy sources for the mass-loss, with 
the bottom (planet interior) pushing material up, and the 
top (stellar heating) pulls the material further away from 
the planet. 


2.2. Race between Cooling and Dispersal 

The importance of our new mechanism depends criti¬ 
cally on the size of the planet when it was first exposed to 
direct stellar irradiation. A highly inflated planet, with 
the photosphere at its Bondi radius, will experience a 
remarkable mass-loss rate of: 


M p ~ 47ri?, B p sur fC s — 47 t Rq 


,2 -Tsurf 


1 x 10 2 yr 


-l 


M„ 


10M fl 


2.0 km s 


0.1 cm 2 g 


-l 


-l 


( 4 ) 


where we have taken the photosphere pressure as P SU rf ~ 
g/n with g being the local gravitational acceleration, k 
the opacity and the chosen sound speed corresponds to a 
separation of ~ 0.1 AU for molecular hydrogen. Clearly, 
such a mass-loss rate is not sustainable and it will bring 
dramatic changes to the planetary atmosphere, an is¬ 
sue we develop later in this paper. For the moment, we 
only consider the on-set of mass-loss. The mass-loss rate 
drops off exponentially with smaller planet sizes, as the 
so nic poin t density depends exponentially on planet sizes 
(Cranmer 2004). So the importance of our new process 
depends on the size of the planet when the disk is re¬ 
moved. 


In turn, how inflated a newly emerged planet is, de¬ 
pends on the race between the cooling/contraction of its 
envelope and the external process of disk dispersal. If the 
cooling time-scale is much shorter than the disk disper¬ 
sal time-scale, then the planet will contract significantly, 
rendering post-partum mass-loss insignificant. 

The relevant cooling time-scale of a gravitational 
bound object with no internal energy sources (e.g. nu¬ 
clear burning) is given by the Kelvin-Helmlroltz time- 
scale (£kh), which for a planet with envelope mass (M env ) 
less than the core mass (M c ) is: 


£kh 


GM c M e] 

R p L 


( 5 ) 


where L is the planet’s cooling luminosity. For an em¬ 
bedded planet, one can assume that any cooling flux that 
reaches the Bondi radius is instantly removed by the ex¬ 
ternal shear flow. This is similar to a planet exposed to 
vacuum. So the actual cooling luminosity only depends 
on the internal structure of the planet, not on its envi¬ 
ronment. 

The near-surface temperature structure determines the 
cooling luminosity. In the deep atmosphere, heat is 
mostly advected outward by convective eddies; near the 
photosphere, this could be accomplished by radiation dif¬ 
fusion. Thanks to stellar heating, the surface temper¬ 
ature of the planet is equilibrated to that of the local 
blackbody (eq. [2| if it is exposed, or to the mid-plane 
disk temperature if it is embedded. This surface searing 
tends to flatten out the surface temperature gradient to 
become more isothermal, reducing the heat escape. This 
situation is similar t o that for an irradiated hot jupiter 
( Guillot e t al. ||l 996[ |Burrows et al.||2000[ |Arras & Bild- 
sten||2UL)5 r 


The thickness of the radiative zone depends on enve¬ 
lope mass, among other things. The planet atmospheres 
are fully adiabatic if their internal entropy is as high as 
the photospheric value- these are the hydr ostatic primor¬ 


dial atmosphe res as calculated by Pollack et al. 


Rafikov (2006). Because of their large scale height, they 


(1996) 


contain the least possible amount of mass. To load more 
mass into these atmospheres, the interiors have to cool to 
a lower entropy than the surface value. The lowering 
of the internal temperature, allied with the surface sear¬ 
ing, then entails an ever thickening isothermal blanket. 
This isothermal layer is effective at blocking the internal 
heat flow. 

As a result, for planets with the same core mass, size 
and surface temperature, the cooling luminosity is the 
largest in ones that have fully convective atmospheres 
(and the lowest envelope masses), while it decreases 
rapidly in ones with an isothermal upper layer and ac¬ 
cordingly larger envelope masses, as is seen in Fig. |T] ob¬ 
tained for planets with a 5rock core with a surface 
temperature of 900 K^] The Kelvin-Helmholtz cooling 
times, for a wide range of models, are presented in Fig. 
Uj where we vary planet photosphere radius (in unit of 
Bondi radius), envelope mass fraction, surface tempera- 


4 The photospheric pressure in these models may differ from 
that at the disk mid-plane, making these models imperfect for em¬ 
bedded planets. However, luminosity from such a model is still well 
characteriz ed as it is largely determined at the radiative-c onvective 
boundary ( [Arras Sz Bildsten|2006~| |Wu Sz Lithwick|2013| >. 
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log 10 (Kelvin-Helmholtz Time-Scale) [years] 








Fig. 2.- The Kelvin-Helmholtz time-scale as a function of planet radius in units of the Bondi radius and envelope mass fraction. The 
solid contours show Kelvin-Helmholtz time-scales of 10 6 and 10 5 years. The left panels show planets with an equilibrium temperature of 
500 K and the right panels show planets with an equilibrium temperature of 900 K. From top to bottom the panels show planets with core 
masses of 3, 5 and 10 M®. 
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ture and core mass. These values are calculate d using 


the MESA stellar and planetary evolution code (Paxton 


et al. 2011 2013), using the setup described in §4. All 


else being equal, one observes that models with a higher 
envelope mass cool much slower. In addition, higher stel¬ 
lar insolation slows down the cooling contraction. Lastly, 
cooling is reduced in models that are less extended in 
size. All of these trends can be explained by the size of 
the isothermal blanket. 

We now turn to discuss the time-scale of disk disper¬ 
sal. There are two time-scales of relevance. The first is 
the time-scale a fully grown planet remains embedded, 
or the remaining disk lifetime after the planet is formed. 
Since we do not know when planets are formed, we adopt 


~ 3- 

10 Myrs (e.g. 

Haisch et al. 

2001 

Hernandez et al. 

2007 

Mamajek et ai.|2009 

). The seconc 

is the time-scale 


As has become clear by now, the two time-scales are 
not necessarily the same. Disk dispersal - particularly in 


the inner disk - is known observationally (e.g., Kenton 
&; Hartmann||1995l |Ercolano et al.|[2011[ iKoepferl et^ al. 
201 3[) to be at least an order of magnitude quicker than 


the disk’s lifetime, the so called “two time-scale” disk 
evolution. These observations have lead to the devel- 
opment of the photoevaporative-switch dispersal model 
(Clarke et al. 2001 Alexander et al. 2006} Gorti et al. 


2009 " Owen et al. 2010 " 2011) where photoevaporation 


of the outer disk carves a gap around 1 AU (for a solar 
mass star), sabotaging the gas supply to the inner region, 
causing the latter to rapidly drain onto the central star 
on its short local viscous time-scale, which is typically 
10 5 years. Alternatively, if the planet is massive enough 
to open a gap, the local disk clearing ti me-sca le can b e 
as short as a few orbital periods (e.g. Crida et al.|2006). 

Having established the time-scales for cooling contrac- 
tion and disk dispersal, we can now assess the size of 
newly emerged planets. Fig. [2] shows that a 5M® 
planet can remain at a size R p ~ Rb well after disk re¬ 
moval, if its envelope mass fraction is more than ~ 10%. 
This threshold reduces to ~ 5% for a 3M® planet, while 
a 10planet can not remain inflated at this size for 
any relevant envelope mass (this is not unexpected as 
10 M® cores a re those typically required for giant planet 
formation e.g. |Rafikov|2006| |Piso et al.|2014| in the outer 
di sk). The last obs ervation is consistent with the finding 
of Lee et al. (2014), where they argue that the envelopes 
of these planets cool and contract so rapidly, they should 
undergo runaway gas accretion even while embedded. 

Further, to this discussion, one may argue that an em¬ 
bedded planetary core continuously accretes more enve¬ 
lope until the envelope cooling time becomes comparable 
to the disk lifetime. Adopt a value of 1 Myrs, this cor¬ 
responds to an envelope mass of ~ 7% for a 5M® planet 
at 0.1 AU. 

In these discussions, we have overlooked an important 
detail. We have assumed that the only planet luminosity 
comes from its cooling contraction. In reality, there is 
the possibility that the planet can accrete planetesimals 
(solids) from the surrounding disk. If this is substan¬ 
tial, it could maintain the high entropy of the gas enve¬ 
lope and keep the planet inflated, well past its Kelvin- 
Helmholtz time-scale. As is shown in Fig, fl} reasonable 
estimates for the accretion luminosity (e.g. Rafikov 2006) 


easily overwhelm the cooling luminosity and may keep 
the planets inflated longer than that estimated in Fig. [2] 
In summary, we argue that newly emerged low-mass 
planets (M < 10M®), fully assembled within the inner 
gaseous disk and have non-negligible envelope masses 
(~ 5 — 30%), will remain inflated ( R p ~ Rb) when the 
disk disappears in a rapid timescale of ~ 10 5 yrs, and are 
thus susceptible to strong winds. sThis stage of mass-loss 



it schematically in §3 and numerically in §4. . 

3. THE “BOIL-OFF” AND ITS CONSEQUENCES 

We now consider the rate of mass-loss and the cooling 
contraction for a low-mass planet, irradiated by a strong 
stellar flux after it emerges from the nascent gas disk. 


3.1. Mass-loss rates 

A highly inflated young planet, exposed to the stellar 
radiation, is vulnerable to mass-loss. If the wind remains 
isothermal, as we have assumed, then the appropriate 
velocity and density profile is given by the Parker wind 
model (Parker| [i958 ). 

Eq. (TTJ), generalized for any planet size R p with R p < 
Rb, becomes: 


M — 4:7rR v p sur fU sur f — 4:7rRpAAp 


Usnrf \ 

C s ) 


4tt GM r 


KC S 


M 


P > 


( 6 ) 


where A4 p is the Mach number of the flow at the planet’s 
photosphere and is only a function of the ratio R p /Rb, 
given by: 


M p = yJ-W 0 [-f (.Rp/Rs )] 


R P 

Rb 


exp - 


2 Rb \ 
R p ) 


(7) 

when R p <C Rb , (8) 


where Wo is a real branch of the Lambert function (Cran- 
mer||2004 1 and 


f(x) = x 4 exp 3 — 


(9) 


We verify that the energy required to heat up the wind, 
given the mass-loss rate (eq. [6]), is subordinate to 
the stellar irradiationfor an equilibrium temperature of 
1000 K, as is shown in Figure [3] for a 10 and 5 earth 
mass planet. So energetically, the wind can be heated to 
isothermal by the star. This energy argument, however, 
breaks down at an equilibrium temperature of ~ 250 K 
for a 5 earth mass planet and ~ 200 K for a 10 earth mass 
planet. Therefore, at large separations, > 1 AU around a 
Sun-like star, the flow cannot remain isothermal, leading 
to lower mass-loss rates. 

Eq. ([6} can be used to calculate a mass-loss time-scale 
as: 


£ml = 


Me ny 

M 


x Ad -1 
4 nG envMp 
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Fig. 3.— Figure showing the ratio of mechanical luminosity in 
a Parker wind, compared to the recived energy from radiation as 
a function of planet radius. The solid line shows a 10 M 0 planet 
and the dashed line shows a 5 planet both with an equilibrium 
temperature of 1000 K. The mechanical luminosity is considerably 
smaller than the energy from irradiation. 


s 10 3 yre X m M ~ 1 ( ) 


^ 0.1 cm 2 g 


where X en v is the envelope mass fraction. In Fig. [4]we il¬ 
lustrate how the Mach number at photosphere, as well as 
the mass-loss time-scale, depend on the ratio of R p /Rb, 
by solving the exact Lambert function. But here, we pro¬ 
vide a schematic argument for the core result. Using the 
expression for the Mach number when Rp < fig (eq. 
[ 8 ]), we can recast the above equation as: 


R P 

Rb 

( T * 


exp 


eq 


\886K 


(2 R b \ 

V Rp J 

- 1/2 


X, 


-1 


0.1 cm 2 g 


-l 


£ml \ 


10 3 years 

-l 


1 


(id 


So at small R p /Rb , the mass-loss is exponentially sen¬ 
sitive to the value of R p /Rb- Alternatively, one can say 
that the final radius where mass-loss effectively stalls is 
only logarithmically sensitive to the exact model param¬ 
eters, and is always of order R p /Rb ~ 0 . 1 . As can 
be observed in Fig. |4j this is indeed the point where 
the mass-loss time-scale sky-rockets to exceptionally long 
time-scales. 

So, once the confining pressure of the disk is lifted, 
the planet’s atmosphere - although in roughly dynam¬ 
ical equilibrium - will undergo rapid mass-loss. Such a 
mass-loss, if happening at a pace faster than the intrinsic 
cooling contraction of the planet, as is indeed the case 
in Fig. [4j will dominate the thermal evolution of the 
planetary atmosphere. 


3.2. Consequences for thermal evolution 

Just like evaporation from the skin cools us in the sum¬ 
mer heat, the rapid mass-loss from young planets will 
lead to dramatic cooling and contraction. 

In a planet without mass-loss, its internal heat is trans¬ 
ported outward by a combination of convection and ra¬ 
diative diffusion. The heat leakage, and hence the rate of 
cooling contraction, is limited by the bottle-neck in this 
transport process, i.e., the outer radiative zone. This is 
the reason behind the results shown in Fig. [5] Mass-loss 



Fig. 4. — Mass loss time-scale (left Y-axis) and launch Mach 
number (right Y-axis) as a function of planet radius normalised 
in terms of the Bondi radius. Thick red lines show the mass-loss 
time-scale for models with k = 0.01 cm 2 g —1 whereas thin black 
lines show the mass-loss time-scale for models with k = 1 cm 2 g —1 ; 
dashed lines indicate models with a = 0.1 AU and solid lines with 
a = 1.0 AU. The thin blue dotted-dashed line shows the launch 
Mach number. We note the two Y-axes are not linked directly. 

alters this picture. The outflow advects thermal energy 
at a rate: 


Tadv ' 


7 -Mel - i n 2 6 —„„-l 


7-1 


10 ergs 


3 x 10 5 ems -1 



where 7 is the ratio of specific heats. This flux easily 
overwhelms the cooling luminosities of low-mass planets 
( p ig- 0 - As a result, the hot interior of the planet can 
shed its energy at a rate much higher than that permitted 
by the radiative bottle-neck. The planet cools down in a 
hurry. 

This cooling reduces the thermal support for the at¬ 
mosphere and it contracts. The contraction reaches 
R p ~ 0.1 Rb within a few /ml- This effectively shuts 
off the mass-loss and cooling resumes the original pace 


& Wu 

2013 

Lopez & Fortney 

2013 

). Unlike the case 

of no mass-loss, where fKH is a 

ways 

of order the age of 


the planet, we now have a planet which is young in age, 
but “looks” old. Specifically, it has a Kelvin-Helmholtz 
time-scale that is considerably larger than its age. 


3.3. Final envelope mass 


Knowing the initial and final radius of the planet en¬ 
velope, we can estimate the amount of mass-loss using 
an energy argument. 

For the initial state, the binding energy of the envelope 


is: 


Ui = 


GM c M l em 

cuRb 


(13) 


where Aj is an order unity constant and is determined 
by the central concentration in the envelope, a is a pa¬ 
rameter (< 1 ) that represents the initial radius of the 
planetary atmosphere in terms of the Bondi radius. Af¬ 
ter the episode of mass-loss, the envelope has shrunk to 
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a radius of 0.1i?s and the binding energy of the envelope 
is now: 

(14 , 

Kb 

where A e is also an order unity number that describes the 
central concentration of the new state. As the mass-loss 


so that the radiation field and the fluid are in local ther¬ 
mal equilibrium. Additionally, we assume that the radi¬ 
ation transport takes place in the static diffusion limit, 
i.e., v/c <C A p /£ (where X p is the mean free path of the 
photons and £ is the scale length of the atmosphere), a 
condition easily satisfied in our problem. 

4.1. Numerical method 


that most of the loss occurs when the planet was at its 

We use the mesa code 

jPaxton et al. 2011 2013 

to 

initial size, and this requires a binding energy release of: 

solve this radiation-hydroc 

iynamics problem. This code 
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(15) 


Assuming that the potential energy difference between 
the initial and the final state is all used to drive the 
Parker wind, or, Ui — U e = U\ os t, we can solve for the 
final envelope mass as: 


M f =( A * + 1 
env VlOaA e + l 


env 


0.2a- 


'Aff 
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(16) 


where we have taken Ai ~ A e ~ 1. 

In this discussion, we ignore two other energy terms. 
One is the energy gain from the stellar insolation, which 
can be orders of magnitude larger than the planet’s in¬ 
ternal luminosity; and one is the energy loss from con¬ 
vective/radiative transport. The latter, as we have ar¬ 
gued, is subordinate to the advective energy flux due to 
mass loss and can be safely ignored (also see numerical 
evidences in ®. The first assumption, however, takes 
some considerations. 

First, stellar heating is essential to drive the Parker 
wind. It continuously heats the gas above the photo¬ 
sphere, allowing them to reach escape velocity. Without 
it, the atmosphere will happily remain at a hydrostatic 
equilibrium. However, stellar heating is only a minor en¬ 
ergy source for mass-loss, not the driving engine. Stellar 
heating provides the final push to take a gas parcel from 
the photosphere to infinity, and the pressure loss caused 
by the outflow allows gas to rise up from below. The 
increase in its binding energy is provided by the gravi¬ 
tational contraction of the underlying gas. If we include 
the energy gain due to stellar irradiation into our simple 
analysis, C/; rr = e(i?s/a) 2 i*tboib where F* is the stellar 
luminosity, tboii is the length of the boil-off phase and e 
a small pre-factor < 1, we find, 




0.2 a 


(17) 


a ) 


GM C 


This shows that stellar irradiation can increase the mass- 
loss slightly, with the strongest effects for small separa¬ 
tions and low core masses. 

In numerical simulations, we frequently find that the 
final envelope mass is some 10% or less of the initial 
mass (Fig. [5]), indicating that stellar irradiation is a sub¬ 
dominant energy source. 

4. NUMERICAL CALCULATIONS 

We simulate the thermal and dynamical evolution of 
a spherically symmetric planetary envelope, subject to 
the boundary conditions of stellar irradiation and mass- 
loss. The gas is optically thick to radiation (t > 2/3) 


traces the evolution of any mass shell inside the planet, 
as is dictated by the equations of mass conservation, mo¬ 
mentum conservation and energy conservation, subject 
to suitable boundary conditions. This is valid so long 
as the simulated domain contains no shocks and sonic 
points. In other words, we can not follow the Parker wind 
through its escape process. Instead, we have to indepen¬ 
dently specify a mass-loss rate at the outer boundary of 
the planet. Since MESA is essentially a Lagrangian hy¬ 
drodynamics code (namely, following the mass element), 
this is not straightforward to implement and we discuss 
this below. 

4.1.1. Implementation of mass-loss 

To specify the mass-loss rate in MESA, one removes 
a certain amount of mass from the top layers (we take 
the surface to be the photosphere with r = 2/3 to the 
outgoing radiation) at the beginning of each time-step 
and then lets the system r eadjust, as is descr ibed in detail 
in the MESA code paper (Paxton et allj |2011 ). 

In our implementation, we assume that the flow above 
the photosphere follows the isothermal Parker wind solu¬ 
tion with a temperature equal to the equilibrium temper¬ 
ature T eq (Eq. [2])and a mass-loss rate as in eqs. (|6|-(|7|). 

At early times when the planet is highly inflatea (it)/ ~ 
Rb)i we find this rate may formally exceed the so-called 
“energy-limited” rate, the rate at which all stellar flux 
is converted into kinetic energy in the wind. For these 
short periods, we cap the mass-loss rate at 10% of the 
energy-limited rate. This choice makes little difference 
to the evolutionary tracks. 

One also has the freedom in how to implement the 
mass-loss term in detail, either over a few grids near the 
photosphere, or spread it across the surface scale heights. 
We have experimented with this and found that provided 
the loss term is smooth, small and does not remove signif¬ 
icant mass from where the contraction of the atmosphere 
is taking place, it does not matter. 

4.1.2. Implementation of flux boundary condition 

The outer boundary is set at the r = 2/3 surface to 
the outgoing radiation. To include stellar irradiation, we 
adopt the MESA F* — E implementation (Paxton et al. 
2013J. This consists of depositing a total irradiatiorTflux 
of F*/4 uniformly down to a column density E, where we 
choose E = 250 g cm -2 as is appropriate for an opacity 
of n v = 4 x 10~ 3 cm 2 g -1 to the incoming stellar flux, 


the value suggested by Guillot (2010) for sun-like stars. 

The F* — E method of including irradiation in planetary 
evolution has been shown to provi de good agreem ent to 
other boundary conditions in |Paxton et al. (2013), when 
mass-loss is not dynamically important. It differs from 
the usual approach of applying a grey or semi-gray ap - 
proximate solution to the atmosphere (e.g. GuillotpOlO), 
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as used by Owen & Wu (2013). It is necessary in this 
problem because energy balance in the upper layers of 
the atmosphere must be taken into account explicitly, 
including PdV cooling and advection. However, for the 
sake of simplicity, we identify numerical models by their 
respective T eq , not their F* values. 


4.1.3. Initialization 

We initialize our models as hydrostatic planets, with a 
radius of Rb (ct = 1), with a chosen envelope mass and 
core mass. The planet models were created by building a 
slightly larger than required planet (at a higher entropy) 
that was then allowed to cool to the required value. The 
core radius is determined assuming a pure rock compo¬ 
sition that fo llows the mass-radius relation from Fortney| 
et al. (2007). We i nclude radioactive h eating from the 
core as described in Owen & Wu (2013) but it has neg¬ 
ligible impact; moreover, we neglect heating from the 
core’s heat capacity as over the short time-scales of in- 
tere st the core is unable to transfer heat to the envelope 
(see Lee et al. 2014 for a detailed discussion). Specifying 
these parameters uniquely determines the properties of 
the planet’s atmosphere. We then switch on the hydro- 
dynamic terms and mass-loss and let the models adjust 
to a new equilibrium (typically this results in a small re¬ 
duction in radius). We then evolve the model forward 
using a time -step that sati s fies t he convergence criterion 
described in Paxton et al. (2011) as well as limiting the 
relative change in envelope mass to < 10~ 3 each time- 
step. We stop the evolution after 3 Myr, at which time 
mass-loss has practically ceased in all our calculations. 
Given the number of simplifications made in order to 
perform the calculations, we perform a series of tests in 
the Appendix to ensure that our results are robust. 

We consider models with two different initial enve¬ 
lope mass-fractions, 10% and 30%, compatible with those 
found by core-accretion models at small separa tions (e.g. 
Bodenheimer fe Lissauer|2014 Lee et ah]|2014 ) and have 
long enough cooling times to remain bloated even at disk 
dispersal. We note the true mass fraction, at the time 
of disk dispersal, is likely determined by the fluid me¬ 
chanics and the thermodynamics around an embedded 
planet and is beyond the scope of this work. We study 
three core masses of 3, 5 & 10 and two equilibrium 
temperatures of 500 & 900 K. 


4.2. The Results 

As is shown in Fig. [5] the evolution of all models is 
similar: the planet’s atmospheres contract, lose mass and 
cool on a time-scale of order 10 s years. The internal 
contraction provides energy to lift up the outer envelopes. 
The mass-loss drops off exponentially with planet size as 
discussed in Section 3.1, and after a few 10 5 yrs, when 
the planet’s radius reaches ~ 0.1 Rb (Fig. [o]), mass-loss 
ceases entirely and the planet cools down normally. At 
this point, the final envelope mass is ~ 10% of the initial 


one, consistent with the estimate in Eq. (16). 

We focus on one model to gain more insight. Fig. [7] 
presents the evolution of a planet with a 5core, an 
initial mass fraction of 10% and at an equilibrium tem¬ 
perature of 500 K (the red solid line in Fig. [5]). Compared 
to a similar model without an imposed mass-loss, the 
contraction of our model is much more rapid, indicating 



Fig. 5.— Mass evolution for models with a 5 M® core. The red 
lines are for models with an equilibrium temperature of 500 K while 
the black lines are for 900 K. The solid and dashed lines stand for 
models with initial envelope fractions of 10 and 30% respectively. 
All models start with a radius R p = Rb- 



Fig. 6.— Evolution of the same models in Fig. [5] shown in the 
radius-mass plane, where the planet’s radius is scaled to the Bondi 
radius. A radius of 0.1 Rb is shown as the thin dotted line. 


that the energy for the mass-loss indeed obtains from 
internal gravitational contraction. The mass-loss time- 
scale (inverse of the mass-loss rate) starts low but sky¬ 
rockets to astronomically long time-scale as the planet’s 
radius shrinks. The Kelvin-Helmholtz time, defined as 
the time-scale to radiate the envelope’s binding energy 
at the current luminosity^] lies always above the system 
age, again indicating that the binding energy is not pri¬ 
marily lost by radiative/convective transport, but by ad¬ 
vection. By the time mass-loss has stopped (~ 1 Myrs), 
the planet has a Kelvin-Helmholtz time of ~ 50 Myrs. It 
has “aged” prematurely. 

We present snapshots of the fluid velocity and lumi¬ 
nosity profiles for the same model in Fig. [8] Initially 
the planet’s envelope is expanding so fast, it completely 
absorbs the planet’s internal luminosity and converts it 
into PdV work, resulting in negative luminosity in the 
surface layers. This layer reaches much deeper than the 
penetration depth of the stellar irradiation, again prov¬ 
ing that the main source of energy for the mass-loss is 
the planet atmosphere’s internal luminosity. By losing 

5 We pick this to be the maximum luminosity in the envelope, 
rather than the surface luminosity as the latter can sometimes be 
negative (Fig. |8|. 
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Fig. 7.— Radius and time-scale evolution for a 5 Mj core and 
an initial envelope mass fraction of 10% , at an equilibrium tem¬ 
perature of 500 K (the solid red line in Fig. |5|. The top panel 
shows its radius evolution with the dashed curve representing a 
similar model where there is no mass-loss. The bottom two panels 
shows how the Kelvin-Helmholtz time-scale (dot-dashed) and the 
mass-loss time-scale (solid) evolve as a function of time and enve¬ 
lope mass fraction. By the time the mass-loss ceases (~ 1 Myrs), 
the planet looks like one that has cooled for ~ 50 Myrs. 

mass, the planet is then able to cool and contract signifi¬ 
cantly such that it shuts off the mass loss. Within 1 Myr, 
the planet’s luminosity has dropped by a factor of 100, 
and its envelope now radially contracts (as opposed to 
outflow). 

Finally, in Figs.[9]&[T0]we present composite results for 
the envelope mass andthe cooling time after 3 x 10 6 years 
of evolution. We see that models receiving a stronger 
stellar insolation lose more mass - roughly by a factor 
of 2 between the 900 K and the 500 K models. We also 
see that less massive planets have a harder time retain¬ 
ing their envelopes. This difference is expected from the 
expression given in Eqn. (17), where stronger irradiation 
and lower core masses lead to slightly more mass-loss, al¬ 
though we note that since the this effect is small and the 
planets cool significantly, energy input from stellar irra¬ 
diation is a sub-donrinant effect. In some cases, envelopes 
around the lowest mass planets (3 M®) can be blasted 
away so completely, they appear as “naked” cores. Fur¬ 
thermore, by the end of 3 Myrs, all planets have a cooling 
time in the range of 4 x 10 7 — 10 s years, approximately 
~ 20 times longer than one would have predicted if the 
planets cooled purely by gravitational contraction. 


5. DISCUSSION 

We discuss how our new process relates to previous 
works on mass-loss of planetary envelopes, and what our 
work implies for the observed planetary properties. We 
connect this process with the first stage of mass-loss, that 


occurs while the disk is dispersing (Ikoma & Hori 


and with the later stage driven by EU V/X-ray p 
vaporation. 


2012 ) 


lotoe- 


5.1. Relation to Ikoma & Hori (2012) 


Ikoma & Hori (2012) considered the evolution of an 
embedded low-mass planet as the background disk den¬ 
sity and pressure slowly declined. They found that the 
planetary atmosphere was eroded, more extremely so if 
an additional heat source keeps the planet puffy. For 
the latter, they invo ked th e coo l ing lu minosity from the 
solid core. However, Lee et al. (2014) argued they may 
have adopted a cooling luminosity that is too large - 
with an assumed cooling time of ~ 10 5 yrs for a core the 
size of the Earth this may have grossly overestimated 
the transport ability of either heat conduction or mantle 
convection. 


Ikoma & Hori (2012) is concerned with the mass-loss 
during the disk dispersal (with a time-scale ~ 10 5 yrs), 
while we are interested in what happens afterwards. So 
the process they investigated is a natural predecessor to 
the one studied here, and it provides the initial input for 
our model. 

Physically, their mechanism is similar to ours, except 
in details, at least in the case of no additional heat 
source . By reducing externa l pressure confinement grad¬ 
ually, Ikoma & Hori (20121 observed a mass-loss that 
likely occurs through a sub-sonic breeze|j In this work, 
we are concerned with a transonic wind. Hoth flows are 
similar in that they both draw energy from the gravita¬ 
tional contraction (when core cooling is not important). 
Above the photosphere, our wind is propelled by stellar 
heating, while their mass-loss is powered by the pressure 
differential between the atmosphere and the background. 

How important is this early episode of mass-loss? or, 
in other words, when does the breeze stop and the wind 
start? We speculate that the transition may occur well 
before the disk becomes so optically thin that the outflow 
can be heated by the star. The infra-red glow from the 
background disk may be already sufficient to keep the 
outflow roughly isothermal. Moreover, the decay of the 
disk’s confining pressure becomes irrelevant once the ram 
pressure in the outflow is strong enough. This occurs be¬ 
fore the gas mean free-path bec omes as large as the plan¬ 
etary radius, a threshold set by 
the termination of their Simula 


Ikoma & Hori (2012) for 
ions. However, in order 


to model thi s properly, one need to improve upon the 


modelling in Ikoma & Hori 
drodynamics terms explicit 


(2012), by including the hy- 
y in the evolution equations, 


and to include the breeze solution explicitly. Ultimately 
a single simulation that follows the protoplanet all the 
way through disk dispersal and boil-off is required to as¬ 
sess the relative importance of the processes in shaping 
the planet population. 

We note that in some of our numerical calculations 
the mass-loss time-scale is similar to the disk clearing 
time scales and as such our results will be sensitive to 
when this cross over takes place. Specifically, if disc 
pressure drops quicker than the “breeze” solutions can 
transfer mass sufficiently to maintain equilibrium, then 
the planet will launch a Parker wind early; however, if 
the disk pressure drops slowly enough, the “breeze” can 
remove sufficient material to make the wind phase (boil- 
off) less important. We note that the breeze mass-loss 

6 However, this process is not explicitly modelled in their work 
and as a result, their mass-loss rate may not be physically self- 
consistent. 
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Fig. 8. — Snapshots of the velocity and luminosity profiles as a function of radius for a model with a 5 Mgj core, an initial mass fraction 
of 10% and at an equilibrium temperature of 500 K (the model shown in Fig. [TJi. From top to bottom, the following times are shown: 10 2 , 
10 3 , 10 4 , 3 X 10 4 , 6 X 10 4 , 10 s , 3 X 10 5 , 10® & 3 X 10® years. Solid lines represent negative values, while dashed lines represent postive 
values. A negative luminosity represents an inwardly directed radiative flux. The large luminosity spike at the surface arises as we are 
explicitly adding the heating from stellar irradiation and is representative of the re-radiated stellar flux. The change in sign close to the 
planet’s surface is abrubt and unresolved, as a result of our simplified F* — £ method for including stellar irradiation. However, as discussed 
in the appendix this simplified treatment does not impact our calculations. 


rates are always lower than those driven by the Parker 
wind, thus we suspect the “boil-off” phase will dominate 
the total mass-loss. 

5.2. Relation to EUV/X-ray driven evaporation 

More mass-loss can occur at the end of our simula¬ 
tions, after the planet has contracted to a radius of 
R p < 0.1 Rb- The ionizing radiation (EUV/X-ray, as 
opposed to continuum radiation considered here) from 
the star can elevate the temperature in the upper layers 
to ~ 5000 —10 4 K, as opposed to the blackbo dy tempera¬ 
ture considered here (Owen & Jackson|2 012). The corre¬ 
sponding sound speed can approach the escape velocity 
at the photosphere, even for a more compact planet. This 
then drives a photoevaporative outflow. This has been 
shown to be particularly significant for Kepler planets in¬ 
ward of ~ 0.2 AU from the host stars a nd can exp l ain the 
smaller sizes of these clo se-in planets (Lopez & Fortney 
2013 Owen k, Wu|2013 ). The atmosphere of these plan- 
ets can be completely removed, consistent with t he hig h 
bulk densities observed in objects like CoRoT-7b (Hatzes 


et al. |20iTj), Kepler- lOb ( Batalha et al.| 20 TIj), Ke jler-3 6b 


(I Carter et al.|[2 012) 

' jOwen &rWuTf50: 


and others (Weiss &; IVlarcy 


2014). 


12013) demonstrated that this mode of 


mass-loss is most significant in the first ~ 100 Myrs, 
mostly because this is when the host star is more chro- 
mospherically active (and hence stronger ionizing radia¬ 
tion), but also because this is when the planet is more 
extended (with a smaller surface escape velocity). Short 
of knowing the initial radii for the planets, they explored 
models with a range of initial cooling times (from 3 Mys 
to 100 Myrs), where a shorter cooling time corresponds 
to models with a more extended photosphere. They con¬ 
cluded that models with the shortest cooling time lose 
the most mass, with the difference being more si gnifi¬ 
can t for lower mass envelopes (Fig. 4 of Owen fe Wu| 


2013). 

This provides the context for our results in Figs. [9} 
[lOl They quantify the initial conditions for the later 
pn otoevap orat ive flo w, just like results from a study like 


Ikoma & Hori (2012) can do for our work here. The early 
episode of mass-loss, propelled by continuum stellar ra- 
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Fig. 9.— Final envelope mass as a function of core mass. Red 
symbols have an equilibrium temperature of 500 K, whereas black 
are 900 K. Open symbols initially had 30% envelope fraction 
(shown as the thin solid blue line) and the filled symbols initially 
had 10% envelope fraction (shown as the thin dotted blue line). 
Assumed to start at a radius oi Rb, all models lose substantial 
amounts of envelopes, with the effect being more drastic for plan¬ 
ets with a lower core mass, a stronger irradiation and a less massive 
envelope. 
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Fig. 10.— Same as Fig. [9] but for the cooling time at the end of 3 
Myrs. The “boil-off” removes internal energy from these planets, 
allowing them to cool off much faster than their ages indicate. 


diation and powered by gravitational contraction of the 
planet, cools the planet to a smaller radii, affecting the 
outcome of later photoevaporation. As such, we argue 
that a coherent study where one follows the evolution 
of the planet envelope through these different stages, in¬ 
stead of studying them independently, would be needed 
to predict the final envelope mass one observes. 

5.3. Applications to Kepler Planets 

Despite being only one of the three possible steps in 
mass removal, the “boil-off” process discussed here may 
explain some of the observed features in Kepler planets. 

It may help explain the dramatic deficit of planets with 
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Fig. 11.— Same as Fig. [9] but for the final planet radius after 10 
Gyrs of cooling contraction. For planets less massive than 10Mq, 
the final radii tend to cluster around the radius observed for 

a large number of Kepler planets. Here, the dashed line indicates 
the radius of the planetary core. 

sizes above ~ 2.5in the Kepler catalogue. This radius 
corresponds to an envelope mass of ~ 1% for planets that 
have cooled for a few Gyrs. Limited by the transit tech¬ 
nique and the mission duration, this catalogue mostly 
contain planets inward of ~ 0.5 AU (T eq ps 400 K). For 
these planets, Fig. [9]shows that, even if they were started 
with much more massive envelopes (10 — 30%), the “boil- 
off” would have left them with much punier envelopes (of 
order 1%) after a few Myrs. Moreover, the final envelope 
mass, at the same location, is found to scale roughly with 
the core mass. As a more massive core exerts a stronger 
gravity on the envelope and compresses it more, this cor¬ 
relation leads to a further convergence in final planet 
radii, as is shown in Fig. El. . . 

In contrast, planets outside this distance range will 
experience less mass-loss (see Fig. [5]) and they may re¬ 
tain most of their primordial envelopes. In other words, 
we predict that there are relatively more neptunes, with 
envelope mass fractions of tens of percent, outside the 
> 1 AU range. This could be tested by future transit 
missions. 

We can contrast these predictions of the boil- off against 
t hose from EUV/X-ray photoevaporation (Owen & Wu 
2013). The latter effect is mostly limited to separation 
^ 0.2 AU. If the latter acts alone, one would expect an 
abundance of Neptunes outside 0.2 AU, closer than our 
above prediction for the “boil-off”. This difference may 
already be discernible in the Kepler data. However, de¬ 
tailed comparison between theory and observation needs 
to take into account of the planet mass, which is largely 
unknown at the moment, especially for planets outward 
of 0.2 AU. 

Are there Neptunes within 1AU? We find that Kepler 
data show a smattering of neptune-like planets within 
this distance. How could these planets have escaped the 
“boil-off”? We speculate that their cores are more mas¬ 
sive than > 10 M ffi (see Fig. |9| and therefore can hold on 
to their envelopes better. This conjecture is supported 
by the general mass-radius relationship of M oc R, as is 
obtained from TTV and RV mass measurements flWu fe 

; Me 


|Lithwick|2013| |Hadden fe Lithwick|2014| Weiss fc Marcy 















12 


Owen, J. E. & Wu, Y. 


2014|. The fact that larger planets do tend to have 


more massive cores is consistent with the “boil-off” pro¬ 
cess. Conversely speaking, the “boil-off” helps explain 
the observed mass-radius relation. 

However, such an explanation is incompatible with the 
presence of a few Neptune-like pla nets w hich have unusu¬ 


ally low densities (Kepler-E 


(Kepler-79, 


Kepler -51, Masuda|| 2014|) 

Jontof-Hutter et ahp0l4|); (Ivepier-87 |(Jtiret al j[20 11 1) . 

„ , „ -fe 


These masses are typically of order 2 — 5Ai®, within the 
range where we predict the “boil-off” should have car¬ 
ried off most of their envelopes. They may be counter 
evidence for our theory, but we hope their mass deter¬ 
minations (all TTV inferences) can be confirmed by ra¬ 
dial velocity studies. Furthermore, since they do not ap¬ 
pear a dominant population (through the radi us distri bu- 


tion) and they may have a different origin (Lee & Chiang 


20151. 


Finally, Rogers (2014) suggest that all planets with 
radii < l.o R® are predo minantly rocky. A 1. 6 Re rock 
planet is roughly a 4 M® (Fortney et ah|2007 |. We note 
that this process may play a role in shaping this result. 
However, we need to be cautious as one needs to perforin 
a more detailed study that factors into th e separation d is- 
tribution, as normal evaporation (e.g. Owen & Wu 2013) 
is able create this as well. Thus, in order to untangle 
the true origin of the solid to H/He transition - which 
may turn out to be a separation dependant statement - 
further modelling is required. 


6. CONCLUSION 

We now know that there are three stages of envelope 
mass-loss: while the planet is embedded in the disk; soon 
after the disk disperses; and while the planet is contract¬ 
ing, with typical durations of ~ 10 5 Myrs, 10 5 Myrs 
and 10 8 Myrs, respectively. These processes combine 
to sculpt the final radius and mass distributions of the 
close-in low-mass planets. 

We have presented and investigated the second stage, 
dubbed the “boil-off” phase, whereby an irradiated, in¬ 
flated (R p ~ Rb) low-mass planet loses its envelope and 
contracts quickly. Our investigation can be summarized 
as follows: 


2. Due to the characteristics of an isothermal wind, 
this “boil-off” will cease when the planet has con¬ 
tracted to within 0.1 Rb- 

3. Order of magnitude energetic considerations show 
that by this time, the envelope mass has reduced 
to ~ 0.2 /a of its initial mass, where a is the initial 
radius in unit of Rb- 


4. We find that this process leads to extreme cool¬ 
ing of the planet. Within a Myrs, the planet have 
cooled down to a state that is only reached after 
~ 100 Myrs of evolution without mass-loss. 


5. Relatively speaking, mass-loss is more severe for 
planets with lower core masses, and for planets sub¬ 
ject to stronger stellar irradiation. Planets more 
massive than 10M® should be less affected by this 
process, so do planets that have accreted so much 
gas they have undergone run-away gas accretion. 


6 . This process could remove the massive envelopes of 
neptune-like planets that are inward of ~ 0.5AU, 
explaining the observed steep fall-off in planet num¬ 
bers above a size of 2.51?®, provided planet for¬ 
mation results in initial envelope mass fractions of 
order tens of percent, as is suggested by recent sim¬ 


ulations and basic time-scale arg uments, (e.g. Lee 


et al. 2014 Lee & Chiang 2015). A corollary is 
that we expect neptunes to be more common at 
larger separations. 


7. This process may also be responsible for the ob¬ 
served mass-radius relation (M ex R). More mas¬ 
sive planets can retain more of their envelopes, 
leading to larger sizes. However, EUV/X-ray pho¬ 
toevaporation can also lead to the same trend. 
Work is needed to identify the relative importances 
of these processes. 


To predict the final envelope mass for planets, one need 
to model the formation, the disk dispersal, the “boil- 
off” and the subsequent photoevaporation process, a task 
that may now be within reach. 


1. We argue that many low-mass planets should 
have photospheric radii R p ~ Rb by the time 
disk disperses, due to their naturally long Kelvin- 
Helmholtz time-scales, and perhaps aided by the 
small amount of heat input by planetesimal accre¬ 
tion.. As such, they are subject to the “boil-off”. 
The primary energy source for the mass-loss is the 
binding energy of the planet itself, with the contin¬ 
uum stellar irradiation acting as a “catalyst” for 
the mass-loss. 
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APPENDIX 


TESTS OF NUMERICAL APPROACH 

In order to model the “boil-off” phase with the MESA stellar/planetary evolution code we have made a number of 
simplifying assumptions in order to model the radiative transfer and mass-loss. Therefore, it is important to check 
that these approximations are not driving our results and our results are robust to the exact tuning parameters chosen 
in MESA. When including mass-loss from the boil-off process in the calculations we assume that the flow is isothermal 
and can be approximated by a Parker wind which we do not attempt to explicitly model. Instead the outer boundary 
for our models is the photosphere to the outgoing IR radiation and we apply the mass-loss as a sink term at the top of 
the envelope. We can check the outer layers of our model (below the photosphere) are quasi-isothermal and close to 
the equilibrium temperature indicating that our choice of a Parker-wind is an appropriate first step. If, for example 
PdV work of the expanding upper atmosphere (below the photosphere) could not be compensated for by the radiative 
heating then the atmosphere would significantly cool below the equilibrium temperature, being closer to adiabatic, 
then a Parker wind would not be an appropriate model. 

Snapshots of the pressure-temperature profile are shown for an evolving planet in Figure |12[ which show the outer 
layers of the envelope are quasi-isothermal during the evolution close to the equilibrium temperature of 900 K. 

As discussed in Section 4, the mass-loss is included as a sub-step in the evolution of the planet’s envelope where 
mass is removed from the upper layers. As noted above, occasionally at early times the mass-loss rate given by the 
Parker wind occasionally exceeded the “energy-limited” rate and we included an efficiency factor of 10% to correct 
this. In Figure 13 we show the evolution of the envelope mass fraction for the model with a 5 core and initial 
envelope mass fraction of 30% where this efficiency is varied with values of 10% (solid) and 100% (dot-dashed), while 
this obviously effects the evolution at early times, the evolution is identical at late times with the final envelope mass 
fraction are within < 0.3% across all the simulations. 
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Fig. 12.— Pressure temperature profiles for the calculation of a planet with an initial envelope mass fraction of 30%, core mass of 5 Mj 
at an equilibrium temperature of 900 K. They are shown at times of 3.7 X 10 4 (solid), 8.3 X 10 4 (dot-dashed), 1.8 X 10 5 (dashed) and 
4.4 X 10 5 (solid) years. 



Fig. 13.— Evolution of the envelope mass fraction for the model with a 5 Mj core and initial envelope mass fraction of 30% at an 
equilibrium temperature of 900 K. The different lines show cases where the mass-loss rates is capped to different values of the “energy- 
limited” rate at early times: 10% (solid - the standard value used in our work) and 100% (dot-dashed). The evolution at early times is 
different, but at late times the evolution is identical and the final envelope mass fractions found are consistent to within < 0.3%. Therefore, 
we are confident the efficiency choice is not driving our results. 


We can also check that the amount of the atmosphere we are smoothing the mass-loss over does not effect our 
results. In MESA the range of the envelope that material is removed every sub-step is controlled by the min_q_f or_k 
parameters that set the mass fraction of the upper envelope over which material is removed. The standard value is 
that any mass is remov ed f rom the top (by mass fraction) < 0.5% of the atmosphere, which is used in all our previous 
calculations. In Figure [14] we show the evolution of the envelope mass fraction for the model with a 5 Mg core and 
initial envelope mass fraction of 30% where the envelope mass fraction material is removed from is < 0.5% (solid), 
< 1% (points), < 2% (dashed) is shown. The final envelope mass fractions at the end of the simulation are within 


< 0.03%. 

In our simulations we use a very simple technique to include the irradiation from the central star in the upper layers 
(but below the photosphere) of the planet’s atmosphere. Where we heat a fixed column density of the atmosphere 
at the rate prescribed by stellar heating, the F* — £ approach (Paxton et al. 2013). We cho se to h eat a c olumn of 
250 g cm -2 in all our above calculations, motivated by the opacity of 4 x 10 -3 suggested by (Guillot]|2010). We can 
check how sensitive our results are to this choice of column density by repeating the calculation" tor a model with a 
5 M® core and initial envelope mass fraction of 30%. We choose to heat a column density of 50 & 1250 g cm -2 at the 
surface of the atmosphere, corresponding to a factor 5 change up and down of the atmospheric opacity. The results 
of this experiment are shown in Figure [l5j As expected the results are obviously sensitive to the depth of the heating 
(as it in turn changes the radius of the planet), the effect of the column density heated by the stellar irradiation is a 
secondary effect and all models follow a similar evolution. 

Considering all the tests above we are happy that our model simplifications are not driving our results, which are 
robust. 
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Fig. 14.— Evolution of the envelope mass fraction for the model with a 5 M® core and initial envelope mass fraction of 30% at an 
equilibrium temperature of 900 K. The different lines show cases where the envelope mass fraction from which material can be removed 
as mass-loss is varied: < 0.5% (black-solid), < 1% (blue-points), < 2% (red-dashed), the evolution is identical in all cases where the final 
envelope mass fractions at the end of the simulation are within < 0.03%. Therefore, we are confident the specific range of atmosphere mass 
over which we remove material is not driving our results. 



Fig. 15.— Evolution of the envelope mass fraction for the model with a 5 M® core and initial envelope mass fraction of 30% at an 
equilibrium temperature of 900 K. The lines correspond to different column densities that are heated by the stellar irradiation, the values 
shown are: 50 (dashed), 250 (solid - the standard value used in our calculations) and 1250 (dot-dashed) g cm -2 . 







